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Abstract 

I report on a numerical program, which can be used to calculate any infra-red safe two- 
jet observable in electron-positron annihilation to next-to-next-to-leading order in the strong 
coupling constant a^. The calculation is based on the subtraction method. The result for the 
two-jet cross section is compared to the literature. 



1 Introduction 



The forthcoming LHC experiment will provide a large sample of multi-particle final states. In or- 
der to extract information from this data, precise theoretical calculations are necessary. This im- 
plies to extend perturbative calculations for selected processes from next-to-leading order (NLO) 
to next-to-next-to-leading order (NNLO) in the perturbative expansion in the strong coupling 
constant. Due to a large variety of interesting jet observables it is desirable not to perform this 
calculation for a specific observable, but to set up a computer program, which yields predictions 
for any infra-red safe observable relevant to the process under consideration. Such a task requires 
the calculation of two-loop amplitudes, a method for the cancellation of infrared divergences and 
stable and efficient Monte Carlo techniques. 

The past years have witnessed a tremendous progress in techniques for the computation of 
two-loop integrals [1-10] and in the calculation of two-loop amplitudes [11-18]. In addition, 
several options for the cancellation of infrared divergences have been discussed [19-29]. Among 
those, the subtraction method - well-known from NLO computations [30-34] - and sector de- 
composition [35-37] are the most promising candidates. For the method based on sector de- 
composition first numerical results have become available for the processes e^e~ 2 jets and 
// — * 7Y [38]. Up to now, no numerical NNLO program which is based on the subtraction method 
and which allows arbitrary cuts and measurement functions is available. 

In this paper I consider the process e^e" — > 2 jets and I report on a numerical program, which 
can be used to calculate any infra-red safe observable related to this process to next-to-next-to- 
leading order. The calculation is based on the subtraction method. The purpose of this article is 
to demonstrate the feasability and correctness of the subtraction method at NNLO in a non-trivial 
example. In the set-up of the program nothing is specific to the process e+e 2 jets - except 
the matrix elements and the fact that in electron-positron annihilation infra-red singularities occur 
only in the final state. Therefore this process - apart from being of interest by itself - serves also 
as a test-ground for the subtraction method for other processes like — > 3 jets or pp ^2 jets. 

I will devote particular attention to the way the phase space is sampled in the Monte-Carlo 
integration and describe methods employed to reduce the statistical error of the Monte-Carlo 
integration. 

This paper is organized as follows: In the next section the subtraction method is reviewed. 

Section 3 describes techniques for an efficient generation of the phase space. In section 4 the nu- 
merical results are presented. Finally, the conclusions are contained in section 5. In an appendix, 

1 collected the explicit formulae for the subtraction terms for double unresolved configurations. 

2 The subtraction method 

The master formula to calculate an observable at an collider with no initial-state hadrons (e.g. an 
electron-positron collider) is given by 
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where pi and p2 are the momenta of the initial-state particles, 2K{s) = 2s is the flux factor and 
s = {pi+ PiY is the center-of-mass energy squared. The factors 1/{2J\ + 1) and 1/(2/2+ 1) 
correspond to an averaging over the initial helicities. d^„-2 is the invariant phase space measure 
for {n — 2) final state particles and On^ {pi, ...,Pn) is the observable, evaluated with an n-parton 
configuration. The index _/ indicates that the leading order contribution depends on ypartons. The 
observable has to be infra-red safe, in particular this impUes that in single and double unresolved 
limits we must have 

oi/\p\,...,p'„) for single unresolved limits, 
0^+2(i^i'---'i^n+2) On\p'i,...,p'J for doublc unresolvcd limits. (2) 

is the amplitude with n partons. At NNLO we need the following expansions of the ampli- 
tudes: 

|2_^(0)*_^(0) , /^_^(0)* (1) (1)* (0)\ / (0)* (2) (2)* (0) W 
\^n.^=^^:^^. (3) 

Here denotes an amplitude with n partons and / loops. We rewrite the master formula eq. (1) 
symbolically as 

(OW) = Y^foi^Uon (4) 

n •' 

and the LO, NLO and NNLO contribution as 













fo. 


^«+2 ""«+2 


+ 1 





^ouy^ = / 0% d<^+ 1 olll d<^l + / oi'> d<^. (5) 

The individual contributions on the r.h.s. of eq. (5) to l^o^i^)^^^ and l^o^i^)^^^^ are in general 
infra-red divergent, only the sum is finite. However, these contributions live on different phase 
spaces, which prevents a naive Monte Carlo approach. To render the individual contributions 
finite, one adds and subtracts suitable chosen terms. The NLO contribution is given by 

(o = / (o% do^l - orf'^)) + / (ol^-) doi^^ + oi/^od^''^) . (6) 

The notation On^ o da^'^^ is a reminder, that in general the approximation is a sum of terms 

Ol^^orfl) = i:ol^)rfl) (7) 
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and the mapping used to relate the n + 1 parton configuration to a n parton configuration differs 
in general for each summand. 



In a similar way, the NNLO contribution is written as 

(oO)r- = /(oW<iai»i-oWo<ia<«J-o«o<iaM) 

+ /KV-;.'i + o^i\o<'aS'>'-oi''»<"«i''") 

+ / (oW rf' + oW orf ^' + oi'' odat'-") . (8) 

d<^^n+\ NLO subtraction term for (n+ 1) -parton configurations, da^n'^^ and dUn'^^ are 

generic NNLO subtraction terms. It is convenient to split these terms into 

au„ - «l*(oo)^„ ""(0,l),n' 

dan - ««(l,0),„ + ««(0,l),n' ^^-^ 

such that (^cxIq q| and da^^^'^^ ^ approximate da^^^2 ^^'^ ^^i+i' respectively, da^^''^.^ ^ and ^i?cx|Q Jj ^ 
are approximations to da^^]^^ . 

2.1 The amplitudes and the subtraction terms 

The NNLO correction to e^e^ 2 jets requires the amplitudes for e+e^ qq up to two-loops, 
the amplitudes e^e~ — > qgq up to one-loop, and the Born amplitudes for e+e qggq and 
e^e~ — > qqq'q' . All these amplitudes are known and can be found in the literature [39^6]. I 
will label the momenta of the particles of the amplitudes as follows: 

Sl[^\\,2) = ^^^\quq2,e+,e^), 

For the subtraction terms I follow the approach of ref. [24] and I use spin-averaged antenna 
functions. Compared to dipole subtraction terms, the use of antenna subtraction terms leads to 
fewer subtraction terms and therefore to a faster program. I follow closely the notation of the 
authors in [24] and denote the three-parton antenna functions by 

4\q,g,q), Df{q,g,g), El\q,q\q'), (11) 

depending on which particles form the antenna. A subscript "sc" or "nf" is used to indicate 
contributions sub-leading in colour or proportional to the number of light flavours Nf. The four- 
parton antenna functions are denoted by 

Af{q,g,gA), Bf\q,q',q',q), cf\q,q,q,q). (12) 
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The NLO subtraction term is rather simple and given by 

j2 



daf'^ = ^AO(l,2,3)|^f(l',20 



(13) 



,(0) • 



The amplitude is evaluated with momenta, which are obtained from the original momenta 
as follows [47] : If {pi,Pj,Pk) is a set of momenta corresponding to an antenna, such that particle 
j is emitted by the antenna formed by particles i and k, then the mapped momenta are given by 



PI 

PK 



{i + p)sijk-2rsjk {\-p)sijk-2rsij 

-Pi + rpj + — N Pk, 



'2-{sijk-Sjk) 
{\-p)sijk-2{ \ - r)sjk 

2{Sijk-Sjk) 



2{Sijk ■ 



^ , {l + p)sijk-2{l-r)sij 



(14) 



where 



r = 



Sjk 



Sij + Sjk ' 



p= /l+4r(l-r) 



Sijk^ik 



For the NNLO subtraction terms, one first needs the NLO subtraction terms for e~^e 
These are given by 



(15) 
3 jets. 



da 



(0,1) 

^^qggq 



\ 1^ [Z)0(l,2,3) +^0(1,3,2) +Z)«(4,2,3) +Z)^(4,3,2)] 



2^[A§(1,2,4)+A§(1,3,4)] 



(0) 



da 



(0,1) 
3,qqq'q' 



1 Nf-l \ N 



4 + 



-[£«(1,3,4) + E?(2,4,3)+£30(3,1,2) + £0(4,2,1)] 



[£0(1,3,2) +^0(2, 4, l)+£3«(3, 1,4) +£0(4,2,3)] jo 



(0) 



• (16) 



,(0) 



Again, the amplitude ^3 is evaluated with momenta obtained from eq. (14), which depend now 
on the partons, which form the antenna. Integration over an one-parton phase space yields 

y ««3 - J [^^3,qggq^^^3,qqq'q') ~ 
1 1 



(0) 



(17) 



The iterated subtraction terms for double unresolved contributions read 



da 



(0,2) ^1/^ 

(0,i),2,^g^9 2 I 2 



[d0(1,2,3)+d0(1,3,2)+d0(4,2,3)+D§(4,3,2)] 
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IN 



da 



IN 

(0,2) 

{Q,\),2,qqq'q' 

+i^[£?(l,3,2)+E?(2,4,l)+E?(3,l,4) + £0(4,2,3)] 



1 , Nf-l\N 
4 2 2 



[£0(1,3,4) +£0(2,4,3) +£0(3, 1,2) +£0(4,2,1)] 



4 2 
N^-l 
2N 



aO(1',2',3') 4^\i",2") 



(18) 



Here we first use the momentum mapping in eq. (14) to relate the four final-state partons to a 
three parton configuration (l',2',3'). The momentum mapping is then used a second time to 
obtain the configuration (1", 2"). The subtraction terms for double unresolved contributions read 



^<oU..^ = i|iV^M(l,2,3,4)+AO(l,3,2,4)] 



1 (NN^-l rn 

J \ \ ^ 



da 



\ N^-\ 
'in 2N 

(0,2) 

Cf 
2 



[a0^,(1,2,3,4)+a0J1,3,2,4)] o 



(0) 



T + K(2,4,3,l)+50(4,2,1,3)) + ^(50(2,4,1,3)+50(4,2,3,1)) 



-^[C0(2,4,3,l)+C0(4,2,l,3)+C0(2,4,l,3)+C0(4,2,3,l)]|o 



(0) 



(19) 



The amplitude Sl^"^ is again evaluated with a two-parton final-state configuration, obtained 
through iteration of eq. (14). For the cyclic order {i,j,k,l) one first compares Sij to Ski. If Sjj is 

smaller, one first combines Pi^Pj.Pk into p'j^p'i^ and then in a second step p'j^p'j^^pi into p'i\p'[. 
The case sj^i < Sij is analogous. Finally we need the subtraction terms for one-loop amplitudes 
with one unresolved parton, They are given by 

^'-Ia0(1,2,3) 



da 



(1,1) _ 

(1,0) ,2 



2N 
N^-l 



(1) 



da 



(1,1) 

(0,1), 2 



4 

2 

N^-l 
2N 



Ai(l,2,3) + ^Ai,„/l,2,3)-^Ai,,,(l,2,3) 



(1) 



[^',{sn) +Nf'E^{sn) + ^^{S23) +%£3°(^23)] - ^-^3 (^o) 



a0(1,2,3) 



(0) 



(20) 



The cancellation of explicit poles of the dimensional regularization parameter 8 occurs individu- 
ally in the following combintations: 



Ja^^i + / da 



,(0,1) 



O (80) , 



"^(1,0) ^"^(0,1) 



O (£«) , 

doi') + J da(i'i) + J Ja(0'2) = o (8°) . 



(21) 



1 2 

Furthermore we have for the process — > 2 jets the additional relation 



(1,1) 

(0,1)- 



(22) 



The correct forms of the unintegrated four-parton antenna functions are listed in the appendix. 
The remaining antenna functions can be found in the literature [24]. 



3 The phase space 

It is a well-known fact, that in the coUinear limit spin correlations remain. For example, the 
spin-dependent splitting functions for g — > and g ^ qq read 



p(0,l) _ 


2 

'ij . 


p(0,l) _ 

g-*gq 


2 



2^^2(1_ 



■4(1-8)Z(1-Z)^ 



where the coUinear limit is parameterized as 

Pi = zp + k 



(23) 



z Ipn 

Pi = i^-z)p-k±_ 



(24) 



l—z 2pn 

Here n is a massless four-vector and the transverse component satisfies 2pkj^ — 2nkj_ = 
The coUinear limits occurs for k^ 0. The term 



1 k I 

^ij K 



(25) 



± 



is proportional to the spin correlation. In four dimensions the spin-averaged splitting functions 
are obtained by integrating over the azimuthal angle cp of around p. By using spin-averaged 
antenna functions, the subtraction terms have not the same point-wise singular behaviour as 
the matrix elements, which is required for local subtraction terms. Instead, cancellations of 
singularities occurs only after an integration over the azimuthal angle over all coUinear splittings 
of the matrix elements. For n final-state particles, this is a one-dimensional integration in the 
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{3n — 4) -dimensional phase space. It can be shown, that in the single coUinear limit, the spin 
correlation depends on the azimuthal angle (p as 

1 

V^-^v ^ Co + C2Cos(2(p + a). (26) 

Sij 

One can therefore perform the average with two points, where the azimuthal angle takes the 
values 

tp, 9 + ^, (27) 

while all other coordinates remain fixed. 



In detail this is done as follows: We partition the phase space into different channels. Within 
one channel, the phase space is generated iteratively according to 

d^n+l = d^„d^DipoleiJ,k (28) 

For each channel we require that the product Sijsjk is the smallest among all considered channels 
and that Sjj < sjk. Therefore it follows that with channel {i,j,k) also channel {k,j,i) has to be 
included into the partioning of the phase space. For the dipole phase space measure we have 

1 1 2k 

d^dipoie = ^l^yi^ -y) I d(p. (29) 



We can therefore generate the {n+ l)-parton configuration from the n-parton configuration by 
using three random numbers wi, M2, M3 and by setting 

y — Ui, Z — U2 (p = 27lM3. (30) 

This defines the invariants as 



Sij = ysijk, 

Sik = z{l-y)sijk, 

sjk = il-z)il-y)sijk. (31) 

From these invariants and the value of cp we can reconstruct the four-momenta of the (n + 1)- 
parton configuration [48]. The additional phase space weight due to the insertion of the (n + l)-th 
particle is 

. = ^(1-.). (32) 

We have therefore a parametrization of the phase space, such that for every coUinear limit the 
azimuthal average can be easily performed, while keeping all other coordinates fixed. It is clear 
that this procedure can be iterated for multiple coUinear emissions. 
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4 Numerical results 



As the nominal choice of input parameters I use N = 3 colors and Nf = 5 massless quarks. I take 
the electromagnetic coupling to be a(mz) = 1/127.9 and the strong coupling to be (Xs{mz) = 
0.118. The numerical values of the Z^-mass and width are mz = 91.187 GeV and Tz = 2.490 
GeV. For the weak mixing angle I use sin^Giy = 0.230. I take the center of mass energy to be 
\/^ = niz and I set the renormalization scale equal to jj^ = Q^. 

As observable I consider the two-jet cross section. The jets are defined according to the 
Durham jet algorithm with y — 0.01. The recombination prescription is given by the E-scheme. 
The two-jet cross section has the perturbative expansion 



(a)(2->0 = (a)(0) ('i + ^fi(2-^0+('^)'c(2-^0y (33) 



In \2kJ 
(a) (0) is the total hadronic cross section at leading order: 

(o)(0) = 40807.4 pb. (34) 
The NLO and NNLO coefficients 5(2-^0 and C^^-J^') have the values 

Bi^-M) ^ -13.674 ±0.004, 

(j{2-jet) ^ -231.6 ±0.3. (35) 

The new result is the value of the NNLO coefficient C^^"]^*), which is obtained directly with the 
methods discussed in this paper. The correctness of this result can be verified with the help of the 
known results for the total hadronic cross section at NNLO, the three-jet cross section at NLO 
and the four-jet cross section at LO. We have the perturbative expansions 

{o)^^ot) ^ (a)(0) + ^5M+(^^j'c(^oO^, 
(a)(3->0 = (a)(0) (^^5(3-^-^^)+ (II) 'c(3->0^. 

(a)(4-^-^^) = (a)(0) (||)'c(4-^-^^). (36) 

The perturbative calculation of the inclusive hadronic cross section (p)^'"*^ is actually known 
to 0{al) [49,50], although we need here only the coefficients up to o(a^). They are given 
by [51-53]: 



^{tot) ^ 2, 

Citot) ^ ^'-1 r/243 ... A,. 3 



8A^ 



4 



= 5.64. (37) 
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The remaining coefficients are obtained by a simple LO calculation (B^^ C^^ -i^'^) or a NLO 
calculation (C^^~j^''>). The values of the coefficients are: 

gi^-jet) ^ 15.679 ±0.004, 
c{3-jet) ^ 153.2 ±0.4, 

^(4-^0 = 84.39 ±0.05. (38) 

These numbers have been obtained with the numerical program presented in this paper. The 
correctness has been checked against already existing programs [48,54]. Since at order ^(a^) 
any event is either classified as a two-jet or three-jet event, and since at order O (a^) any event is 
either classified as a two-, three- or four-jet event, we must have 

(-{tot) _ f-,{2-jet) _^f-,{3-jet) _^f-,{4-jet)^ ^^9) 

We find 

B{2-jet)_^g{3-jet) ^ 2.005 ±0.006, 
Ci2-jet)_^(.{3-jet)_^(.i4-jet) ^ 6.0 ±0.5. (40) 

These numbers agree nicely with the values given in eq. (37) for the total hadronic cross section. 

5 Conclusions 

In this paper 1 reported on a numerical program for two-jet observables in electron-positron 
annihilation at next-to-next-to-leading order. To cancel the infra-red divergences the subtraction 
method with antenna functions is used. The correctness of the numerical program is verified by 
comparing the results for the two-jet cross section to values, which can be obtained indirectly 
by subtracting from the known result of the total hadronic cross section at the next-to-leading 
order result of the three-jet cross section and the leading order result of the four-jet cross section. 
The numerical program is set up such that - apart from the specific matrix elements - nothing is 
specific to the process e~^e~ — > 2 jets. I have therefore confidence, that the subtraction method 
can be extended to other processes like e~^e~ 3 jets or pp ^2 jets. 
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A Antenna functions 

Here 1 list the correct unintegrated four-parton antenna functions 

(?,g,g,^), Afli^,g,g,q), Bf\q,q',q',q), cf\q,q,M)- (41) 
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The formulae printed in ref. [24] differ significantly from these results. However, the formulae 
presented here agree analytically with a FORM-file obtained from one of the authors of ref. [24] . 
For the integrated antenna functions 1 found agreement with the results published in ref. [24]. 
Each antenna is given as a product of a prefactor fp , a symmetry factor S , a colour factor C and 
a kinematical factor ^ . The prefactor is given for the unintegrated antenna functions by 

'P = {SnUsf. (42) 

The symmetry and colour factors are all equal to one, except for the antenna C4, which has a 
symmetry factor of 1 /2. The kinematical factors are given by 



aO(1, 2,3,4) 



1 ^ 4512 + 4534-851234 ^ 2 ^ 4 ^ 2 ^ 2 2512^ + 2534^ 



534 5234 ■^23 ■^HS ■^234 ■^12 ^^34 ■^12 ■^234 ■^23>^123 ^^23 ^^234 ^^23 ■^123 ^^234 ■^1234 

-25123 - 2 5234 + 251234 ^ Sm^ + S34^ ^ 3523 -35123 ^ 3534 -35234 ^ 4512534 



512 523 534 512 523 5234 51234 5125234 51234 523 5123 51234 523^5123 5234 

512 + 5234 ^ -2 523 + 3 534 - 451234 ^ 5i23 +5234 ^ 2 523 51234 + 523 + 251234 



523512353451234 5i2 5123 5234 512 52353451234 5125123 5345234 

4 6 6 6 1 3 3 2 



■ + + + y+ y+ , 

5123 534 5123 5234 5123 51234 5234 51234 5125123 5i23 5234 523 

-3512 + 35123 523^+534^ ^ 3523 -35234 _ -3512 + 2 523 + 451234 



523 5234 51234 5i2 5i23 5234 51234 5i23 534 51234 5123 534 5234 

2523 -5234 + 251234 2523 - 5i23 + 251234 5i2^+523^ _^ -25i2 - 2 5234 + 251234 



512 5123 534 512 534 5234 5i23 534 5234 51234 523 5123 534 

fl2" 



3512 + 6523 + 3534 , 523 , ^534 , 45i2 4 534 , 2534^ , 251"^ 



5123523451234 5345234^ 5235234^ 5235123^ 523^5234 523^5234^ 523^5123^ 

4512 25123 + 2534-251234 ^ 523 ^^2) 
523^5123 512 5235234 512 5123^ 

The subleading-colour antenna function A4 ^^(1,2, 3,4) is split into two parts 

a0,,(1,2,3,4) = Ia0,,(l,2,3,4) + ia0,,(l,3,2,4), (44) 

such that ^,^.(1,2,3,4) corresponds to the cyclic ordering (1,2,3,4), while 04 ^^(1,3, 2, 4) cor- 
responds to the cyclic ordering ( 1 , 3 , 2, 4) . ( 1 , 2, 3 , 4) is given by 



<.c(l,2,3,4)| 



512 + 513 - 2514 + 524+534 _^ -5l2 + 523 - 5l3 - 5l4 _ 513 524 _^ 1 



523^51234 534 5234 51234 523 5 123 534 5 1234 523 51234 

8512514 + 3 512 524 + 512 534 + 8 514 513 + 513 524 + 3513 534 

2 5235123^51234 
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2 523 5123 523451234 

3514513 - 3514523 + 4514^-513523+513^ + 523^ ^ 523-2513-2514-524 



512 534 5234 51234 512 5234 51234 

3512524 + 512 534 + 513 524 + 3513 534 + 8514524 + 8514534 _^ -25i2 - 2 5i3 - 8 5i4 + 524 + 534 



2523 5234^51234 2523 5234 51234 

513 524 ^ 523 - 513 - 2514 - 2524 ^ 2514^ 



523^-53451234 512353451234 5123534523451234 (513 +534) 

2514^ _l_ 2514^ 



51251234 (512+513) (512 + 524) (513 + 534) 51253451234 (512 + 524) (513 + 534) 
2514^ 2514 1 



■ + ■ 



53451234 (512+524) (513+534) (524 + 534) 51251234(512+513) 53451234 
2514^ 2514 1 



+ 



512 5123 5234 51234 (512+524) 53451234 (524 + 534) 51251234 

-2512513 + 4512514 + 512524 - 512 + 4514513 + 513 534 - 513 , 523 (514 + 512 + 513) 



+ ■ 



523^5123 51234 534 5234^51234 

513 524 ^ 513 524 ^ 4514 523 + 524 523 - 523^ + 514 524 + 514 534 

512 523 ■^234 51234 5i2 523^51234 5i2 5i23 5234 51234 

3 5 14 524 - 3 5 14 523 + 4 5 14^ - 524 523 + 524^ + 523 ^ 



512 5123 ■^34 51234 

-512524 - 513534 - 4514524 - 4514534 + 2524534 + 524^ + 534^ , 35i2 + 3 5i3 + 6514 



+ ■ 



523^523451234 2 5234^51234 

6514 + 3524 + 3534 -512 514 - 513523 - 4514523 + 523^ - 514513 513524 (512+513) 



25123^51234 5123 534 5234 51234 523^5123 534 51234 

512 -523 +513 +4514+524+534 _^ (5l2 +5l3) (25i2 5i4 + 5i2 524 + 25i45i3 + 5i3 534) 



5123 5234 51234 523^5i23^5i234 
6514 + 3513 + 3524 513524(524 + 534) , 25i4 (512 + 513) (524 + 534) 



+ 



512 534 51234 5i2 523^5234 51234 523^5123 5234 51234 

514 (512 + 523 + 2514) ^ (524 + 534) (513 534 + 2514534 + 512524 + 2514524) 



534 5234 51234 (513 + 534) 523^5234 51234 

-512-513 + 8514 + 2524 + 2534 ^ 514 (512 + 2514) ^ 514 (534 + 523 + 2514) 



2 523512351234 53451234 (513 + 534) (524 + 534) 512 512351234 (512 + 524) 

514 (513 + 523 +2514) ^ 5i4 (523 + 524 + 25i4) ^ 514(2514+524) 



512523451234 (512+524) 5i235345i234 (513 +534) 5i25345i234 (513 +534) 

514(513 + 2514) ^ 514(513+2514) ^ 514 ( 2514 + 524) 



53451234 (512+524) (524 + 534) 512 534 51234 (512 + 524) 5i2 5i234 (5l2 + 5l3) (5l3 + 534) 

514(534 + 2514) (523 --^u) (523^-2514 523 + 2514^) _^ 523 (5l4 + 524 + 534) 



51251234 (512+513) (512 + 524) 512 5123 534 5234 51234 5i2 5i23^5i234 
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-523 + 514 +524 +534 



512345123 512 

Finally the antenna functions with four quarks: 

50(1,2,3,4)1^ = 

-2 ^14 S24 534 - 512 524 534 - 5l3 524 ^34 + 5l3 524^ + 5l2 534^ 
523^5234^51234 

[2*13 514 524 + 2 512 514 534 - 512 513 524 "512 524 534 "513 524 534 + 513^524 



523^5234512345123 



+513 524^ - 512 513 534 + 512 534^ + 512^534] 

-2 512 513 514 - 512 513 524 + 513^524 - 5l2 5i3 534 + 512^534 ^ 5i4 524 + 514 534 + 5i2 524 + 513 534 

523^51234 5123^ 523 5234^51234 

514 (524 + 512 + 513 + 534 + 2514 ) 5l2 5l4 + 5i4 5i3 + 5i2 524 + 5l3 534 2^14 

H 1 2 1 ' (^"'^ 

523 523451234 5123 523 5i234 5l23 523451234 5123 

C0(l,2,3,4)|^ = 

512514513 ^ 512 (514 -513) _^ -512 514 - 514513 - 514 523 " 514^+512513 +512534 



524 5124 523 51234 5123 2 524 5124 523 51234 2 524 5124 5234 51234 

512514 512(514-513) 512 _^ 534 (514 + 512+513) 



524 5124 51234 5123 2524523 512345123 52452351234 524 5234^51234 

512514 +514513 +514523 + 514^ -512513 + 512534 _^ 5i2 



2524 5234 51234 5123 524 5234 51234 

2 

512 514 -'V14 513 + A-12.V24 -513 A'24 -512 A'13 - 513 512^-13 _^ 5i2 



25124 523 523451234 5124 523512345123 2 5124 52351234 

514 + 512 - 513 512 524 (514 + 512 + 513) 



25124 5234 51234 5i24 51234 5123 523 5234^51234 

-512514 + 514513 + 512524 + 513524 + 512513 + 513^ , 5i4 + 25i2+5i3 , 5i2 



+ - + ■ 



2523 5234 51234 5123 523 5234 51234 2523 51234 5123 

514 + 512+513 ^ 514-512-513 ^^^^ 



5234 51234 2 5234512345123 
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